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Abstract 

The 2D off-critical q-state Potts model with boundaries was studied as a factorizable 
relativistic scattering theory. The scattering S-matrices for particles reflecting off the 
boundaries were obtained for the cases of "fixed" and "free" boundary conditions. In the 
Ising limit, the computed results agreed with recent work [5]. 
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Introduction 

In the past decade, many advances were made in 2 dimensional statistical mechanic 
systems by applying the idea of conformal invariance[l]. Their critical points are well 
described by conformal field theory (CFT) and their universality can be classified by a 
Virasoro central charge c[2}. When such a model is perturbed off their critical point, the 
conformal symmetry is broken and the theory develops finite correlation lengths. However, 
for certain perturbations, residue symmetries survives in the form of an infinite set of 
commuting integrals of motion and renders the theory integrable[3]. Such "perturbed 
conformal field theories" can often be described by a relativistic scattering theory of massive 
particles where the S-matrix is factorizable. All physical information about the field theory 
can then be obtained from the S-matrix by constructing correlation functions using the 
form- factors method [4]. 

Recently, statistical systems with boundaries were studied using the above method [5]. 
It was found that one can choose certain boundary conditions which perserves the integra- 
bility of the bulk theory. Such integrable boundary conditions can be represented by the 
boundary S-matrix describing the scattering of particles from the boundary. In this work, 
such boundary S-matrices are obtained for the q-state Potts model (for < q < 4) with 
free and fixed boundary conditions. 

Integrable q-states Potts Model 

In the lattice q-Potts model [6], the spins a(x) at the sites of the lattice are allowed to 
be in one of the q possible states (1, 2, q). The partition function has the form 

Z = E II (l + KSaWrtv)), (1) 

{a(x)} (x,y) = <nn> 
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which is invariant under the permutation group S q . The phase transition point of this 
model occurs at 

K = K c = y/q. (2) 

For < q < 4, this is a second order phase transition, and its critical point is described by 
the CFT with central charge [7] 

6 , . . 7T p — 1 

c=l ; where Jq = 2sin{ ). (3) 

p(p+l) vy y 2p + l> v ; 

The energy density e(x) of the Potts model then corresponds to the degenerate primary 
field $(2,1) with dimension 

This field is also a relevent operator in the theory, and preturbation of the CFT action by 
this field leads to an off-critical theory with action 

Aq, T = A CFT(c) + t J e(x)d 2 x (5) 



where 



K c -K 

T = -IcT ' (6) 



L c 



and A C ft(c) is the "action" of the CFT with central charge c. This action describes the 
scaling domain of the q-Potts model, and is shown in [8] to be an integrable field theory (ie 
possess nontrival higher-spin local integrals of motion) . Following [8] , we consider the low 
temperature phase K > K c of this theory. The field theory (5) has q degenerate vacua | a) ; 
a = 1, q, with S q acting by permutations of these vacuum states. Its particle content 
must then contain q(q — 1) "kinks" A a },; a, b = 1, q; a ^ 6, which corresponds to the 
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domain wall separating vacua a and b. The antiparticles can be identified by A ab = A ba . 

i 

The masses of these particles are taken to be equal to M ~ | r 1 1_Ae . 

The scattering of these asymptotic particles is governed by the S-matrix and the 
integrability of (5) implies that this S-matrix is factorizable. Recall that the energy- 
momentum of particles can be parametrized by their rapidity 9, where 

E = MchO, P = Msh6. (7) 

The asymptotic particles states are generated by the "particle creation operator" A ab (6) 
satisfying the quadratic commutation relations 

A ab {0i)A hc (e2) = J2 s at(® 12) A ad (6 2 )A dc (6 1 ), (8) 

where 6*12 = Q\ — 6*2- The S q symmetric two-kink S-matrix elements (Fig.l) were computed 
in [8] to be (for a 7^ b 7^ c 7^ d) 



s») = s (e) = (*-,«» n(jg ); (9.) 

sh(\(6 - 2p)) s fe(A(f - 9)) in h 



s a M - s 1 (9) - sm(Af) sMA(¥ _^ n (-)' 



^) = ^)= ri n(A|) jM A(g-V)) n( ^ ); ^ 



where 



and 



r(i-x)r(i-A + x)r(|A-x)r(|A + x) /x , 
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T(l + 2k\ - x)T(2k\ - x)T(l + (2k - §)A - x)F((2k + f )A - x) 
r(l + (2k + 1)A - x)r((2fc + 1)A - x)T(l + (2k - |)A - x)T((2k + |)A - x) ' 

(11) 



The amplitudes Si(9) possess a "bound-state" pole at 6* = ~p (Fig. 2a) while 82(B) have 
the "cross-channel" pole at 9 = ^ (Fig. 2b). Likewise So(#) have both the above poles 
(Fig. 2c) with residues 

Res 0= 2^ 1 So(9) = Res e= 2KiSi(6) = 
-Res e=if S Q (6) = -Res e=lf S 2 (6) = if(X), (12) 

where 



(1 - e-^-l x ^)(l - e~§ At )(l - e? xt )e- xt dt 
2(l-e-*)(l + e- At ) T 



} (13) 



is the "three-kink coupling" (Fig. 3). 

As is shown in [8], the complete particle spectrum of field theory (5) for 3 < q < 
4 is quite complicated with the appearance of exotic excitation and bound states. For 
simplicity, we will restrict our attention to the range < q < 3. 

Boundary S-Matrix 

It is natural to consider the Potts model in the presence of boundaries with some 
boundary conditions imposed on the boundary spins. In [5], it was shown that certain 
boundary conditions preserves the integrability of the bulk theory (ie an infinite subset of 
the bulk integrals of motion survives with the introduction of the boundaries) . For a rela- 
tivists scattering theory of massive particles, one can associate these integrable boundary 
conditions with certain boundary S-matrix which describe the scattering of particles with 
the boundary. 



To be more precise, consider the model defined on an semi-infinite plane (say x G 
(—00, 0], y G (—00, 00), the y-axis being the booundary). Let us suppose that there exists 
integrable boundary conditions for the Potts model with the modified action 

/oo 
dy$ B (y), (14) 
-00 

where A qjT+ cBC is the action (5) with certain conformal boundary conditions (CBC), and 
®b(v) is some relevant boundary operator [9,5]. One can think of (14) as a perturbation of 
CBC, and the corresponding Fock states can be classified as asymptotic scattering states. 
In particular, the boundary with boundary spins in the state "a" can be associated with a 
stationary impenetrable particle B a of infinite mass at x = 0. Then the asymptotic n-kink 
scattering state can be written as the product 

(e 1 )A (e 2 )...A 

a n -ia n 

(e n -{)A ana {e n )B a , (15) 

where the vacua ai, a 2 , a n satisfy the restrictions a* 7^ o>i+i, and a n 7^ a. 

If the initial "in-state" of scattering is the asymptotic state (15) with 6\ > 6 2 > 
... > 9 n > (ie n kinks moving towards the boundary of state "a"), then in the infinite 
future, this state becomes a superposition of the final "out-states". Integrability of (14) 
constraints "out-states" to have the form 

A blb2 (-e 1 )...A bn _ lbn (-e n - 1 )A bnb (-e n )B b , (ie) 

with hi 7^ bi+i and b n 7^ b. Thus we have the relation 

A ai a 2 (9i)A a2a3 (92)...A an _ ian (9 n -i)A aria (9 n )B a = 
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E - E E ^V.iK^l, ^)^ bl 6 2 (-^l)...^6„_ 1 6„ (-*„-l) A„6(-^n)5 6 , (17) 

which defines the n-kink Boundary S-matrix. When n=l, we have the simple commutation 
relation 

A ba (e)B a = J2 R ba(0)A bc (-0)B c , (18) 

c 

where R^ a {6) are elements of the boundary S-matrix for the reflection of one particle off 
the boundary(Fig.4). It follows from the factorizability of the scattering that R^\ 
can be written as a product of bulk amplitudes S^, and boundary amplitudes R^ a . For 
example, when 2 kinks scatter off the boundary, the amplitude for this scattering can be 
factorized in two equivalent ways(Fig.5), leading to 

E E Rf ba(^)s b c 9 f (0i + e 2 )R e gf {e 2 )sf e {e 2 - e 1 ) 

9 f 

= EE 5 - (fc +fc wo. (19) 

9' f 

which is known as the "boundary Yang-Baxter equation" [10]. 

As is known, the amplitudes Rl a {0) have to satisfy several conditions in addition to 
(19). Firstly we have the "boundary unitarity condition" (Fig. 6) 

Y, R Sa(0)Ric(-0) = 6i, (20) 

c^b 

which is a direct analog to the unitarity condition for the bulk S-matrix. To obtain the 
crossing symmetry condition for the boundary scattering, it is necessary to use the "cross 
amplitude" [5] 

K abc (e)=RU^-9). (21) 
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As shown in [5], this amplitude K abc (6) has to satisfy the so-called "boundary cross- 
unitarity condition" 

K abc {6) = SH{26)K adc {-e), (22) 

d=£a,c 

which is illustrated in Fig. 7. Finally we have the "boundary bootstrap equation" [11] which 
describes the scattering of "bound-state" particles with the boundary. In the bulk theory, 
the kink A ab can appear as a bound-state particle in the two particle scattering (see Fig. 2) 

A ac {6 + j)A cb (9 - j ) - fA ah {6) (23) 

where / is the 3-particle coupling in (13). Applying the algebras (8), (18) and (23) to the 
"in-state" A ac {9 + ^)A cfc (# — ^)B b , we obtained the bootstrap equation(Fig.8) 

KM =EE Rf *V - J^mRtAO + f ) (24) 

for the S q symmetric Potts model. 

Equations (19) through (24) allows the boundary S-matrix elements R c ab (0) to be 
determined up to some CDD factors [5]. For boundary conditions which respect the S q 
symmetry, we can expect R c ab {9) to have a pole at 9 = ^ (Fig. 9a) with residue 

KM^'-M^, (25) 

where g£ is the amplitude for coupling of the particle A cb to the boundary(Fig.lO). Fur- 
thermore, if g%, gl ^ 0, the element R c ab {9) has another pole at 6 = ^ (Fig. 9b) where 

^6(0)^^4; (26) 

this pole is shown in Fig. 10. Of course the presence of the above poles depends on the 
boundary condition, as we shall see when we consider the two simplest cases: "free" 
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and "fixed" boundary conditions. Both cases are conformal boundary conditions (i.e. 
&b(u) = 0) and we conjecture that they also preserves integrability in the off-critical 
Potts model. 

Fixed Boundary Condition 

In this simple boundary condition, the boundary spins are all fixed to one state, say 
"a" . The corresponding boundary S-matrix element 

R a ba (0) = Rf(0) (27) 

satisfies the boundary Yang-Baxter equation (19) automatically. To determine this ampli- 
tude, one appeals to the unitarity condition (20) 

R f (9)R f (-9) = 1; (28) 

and the crossing symmetry condition (22) 

K f (9) = [(q - 2)S 2 (29) + S 3 (29)]K f (-9), (29) 

where 

K f (9)=R f C{-9) (30) 

is the crossing amplitude. Since all boundary states are fixed, we do not expect Rf(9) to 
possess any poles in the physical domain < 9 < ^. The solution to (28) and (29) can 
be f actor ized as 

\9 

R f (9) = F (9)F 1 (-), (31) 
it: 



where F\[X) solves 



sm(27r(A-X)) 
sm(27r(| -X))' 



Fl (X) = B(A - 2X) )J_) X - ^.[[ Fi(X - X); (32) 
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F 1 (-X)F 1 (X) = 1, (33) 
and its minimal solution can be written as 

with 

r[(4/c - 1)A + 2X]r[(4A; - 3)A + 2X + l]T[(4k - 3)A + § + 2X]T[(Ak - 2)A + § - 2X] 
T[4A;A + 2X]r[4(/c - 1)A + 2X + l]r[(4fc - 2)A - J + 2X + l]r[(4fc - 1)A - J - 2X + 1], 

(346) 

up to some CDD factors. 

The factor F (6>) can be obtained from the fixed boundary bootstrap equation (24) 



77T 77T 

Rj{9) = [8,(29) + (q- 3)So(26)]R f (e - -)R f (9 + -), (35) 



which reduces to 



/ ,7T id . . TV %d . . 5"7T i0,^, n *7Ts„/„ ^ . , . 

F O (0) = -ton(- + (- + -)cot(— + -)F O (0 - - )F O (0 + -), (36) 



with simple solution 



F o (0) = -tan(^ + |). (37) 



In the Ising limit (g = 2, A = |), the boundary S-matrix have the form 



R f {0) = itanCj-l), (38) 



as obtained in [5] by field theoretic method. 
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In the other interesting limit q = 3 (A = 1), suppose the boundary are fixed in the 
state "A" , the solution simplifies as 



R fW) = ~ ■ ,1 ie V ( 39 ) 
sm(f - f ) 



Free Boundary Condition 

In contrast to the fixed boundary condition, we have the "free" case where the bound- 
ary spins can be in any one of the q states. The corresponding boundary S-matrix has to 
respect the S q symmetry and the algebra (18) simplifies to 

A ba (6)B a = R 1 (9)A ba (-9)B a + R2(e)A bc (-0)B c , (40) 

c^a,b 

where the amplitudes R\{9) and R2(0) are shown in Fig. 10. 

The boundary Yang-Baxter equation (19) provides three non-trivial equations: 

R 1 S 3 R 2 S 1 + R 2 S 1 R 1 S 1 + (q- 3)R 2 S 1 R 2 S 1 + {q- 3)R 1 S 2 R 2 S 

+(q - 3)R 2 S RiS Q + (q- 3)(g - 4)R 2 S R 2 S Q 

= (q- 2)R 2 S 2 R 1 S 2 + (q- 3)R 1 S R 2 S 2 + (q - 3) 2 R 2 S R 2 S 2 

+R 1 S 1 R 2 S 3 + R 2 S 3 RiS 3 + {q- 3)R 2 S 1 R 2 S 3 ; (41a) 

R\S\R 2 S 2 + (q — o)R 2 S\R 2 S 2 + R 2 S 2 RiS 3 + (q — 3)R 2 SqR 2 S 3 

+(q - 3)R 2 S 2 R!S 2 + R 2 S 3 R t S 2 + (q - 3)R 1 S R 2 S 2 + (q - 3)(q - 4)R 2 S R 2 S 2 

= R\S 2 R 2 S\ + (q- 3)R 1 S 2 R 2 S + (q - 3)R 2 S R 2 S 1 + (q - 3) 2 R 2 S R 2 S Q ; (416) 

RiS 3 R 2 Sq + R\S 2 R 2 S\ + (q — 4)RiS 2 R 2 So + R 2 SiR\So + R 2 SoRiSi 
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+(q - 4)R 2 S R 1 S + {q- 3)R 2 S 1 R 2 S + (q - A)R 2 S R 2 S 1 + (q - 4) 2 R 2 S R 2 S 

= R\S\R 2 S<2 + RiSqR 2 S3 + (q — A)R\SqR 2 S 2 + R 2 SsRiS 2 + R 2 S 2 RiSs 

+{q - 3)R 2 S 2 R 1 S 2 + {q- 3)R 2 S 1 R 2 S 2 + {q- 4)R 2 S R 2 S 3 + [{q - 3) + (q - 4) 2 ]R 2 S Q R 2 S 2 ; 

(41c) 

where the argument in each term has the form Ri{6\)Sj{9i + 9 2 )Rk(9 2 )Si(9 2 — 6>i). 

Equation(41) can be solved for the ratio Ri{9i)/ R 2 {9\) when we take the limit 9 2 — > ^ 
and noting that both Ri(9 2 ) and R 2 (9 2 ) have a simple pole in this limit with the same 
residue 

Res d2 ^R 1 (9 2 ) = Res 02=lf R 2 {9 2 ) = l -g 2 (42) 

for some boundary coupling g. 

The solution to (41) can then be written as 

*.W = fa-3)^ + %P(g) (43a) 
v ' vy sinhA(z7T - 29) y m' v ' 

d m = sin^ sinh2Xe sinhA(f + 2g) Ag 

U2[0) sin f sinhA(z7r-20)sinhA(f -20) 1 ztt J ' K } 

where we use the fact that R 2 (9) has a simple pole at 9 = *jr, which is absent in R\{9). We 

do not expect R\{9) and -R2(#) to have any other poles in the physical domain (0 < 9 < ^ ). 

The normalization factor P(j^) is constrained by the unitarity conditions (20) and (21), 

which reduces to 

P(9)P(-9) = 1, (44a) 



and 

m 9 \9 sinhA(nr + fl) sinhA(^f + 9) p in 9 

{ Y~r~ ^sinhA(f -9) sinhA(^ -9) ( T + 2 
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respectively. The minimal solutions can be written as 
where 

T[(4/c - 1)A + 2X]T[(4k - 3)A + 2X + l]T[4k\ + § - 2X]r[4(£; - 1)A - | - 2X + 1] 
T[4kX + 2X]T[4(k - 1)A + 2X + l]T[(4k - 3)A + § - 2X]T[(4/c - 1)A - J - 2X + 1], 

(456) 

up to some CDD factors. The sign in (45) can be justified by the boundary bootstrap 
equation (24). 

In the Ising limit (A — > |), the boundary S- matrix element simplifies to 

fli(0) = -cot(£ + ^), (46) 

in agreement with [5]. For the q = 3 Potts model (A — > 1), the two scattering amplitudes 
have the form 

i?i(0) = 
R (d) _ sm(fL-f) s in(f-f) 

^ " sin (f L + f ) sin(f + f )- < 47 > 
Finally, a simple computation give us 



/' (3 ~ q) . ,4tt f <lt e"*" 4 ** 

S'(A) = -\ sin(— X)exp{ 



A v 3 ' ^Vo * 2(l-e-*)(l + e- 2At ) 

[(1 _ e -At) (1 _ e 3At- t) + e ^-t (1 _ e -2At+M +t)(1 _ e 3At )]} (4g) 



for the boundary coupling constant. 
Conclusion 
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In this study, the q-Potts model boundary S-matrix for free and fixed boundary con- 
ditions were dervied. It would be interesting to apply the techniques of Thermodynamics 
Bethe Ansatz to study these S-matrices. In particular, one can investigate the renormal- 
ization group flow between these two boundary conditions[5,13]. This work is in progress. 
Finally, we would like to note that it may be possible to use these boundary S-matrices to 
compute the crossing probabilities for the percolation problem (q = 1) in a finite region[12]. 
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Figure Captions 

Fig.l. The scattering processes for the bulk S-matrix elements S (9) : Si(9), S 2 (9) and 
S s (9) defined in (9); with a^b^c^d. 

Fig. 2. Space-time diagrams associated with the pole at 9 = in Si (9) [Fig. (a)] and 
the corresponding cross-channel pole in S 2 (9) [Fig.(b)]. The amplitude S (9) possess both 
poles as shown in [Fig.(c)]; the states a, 6, c, d are all different. 

Fig. 3. Tree-kink vertex associated with the "coupling constant" / 

Fig. 4. The boundary scattering processes described by the amplitude Rl a (9), with 
b 7^ a, a 

Fig. 5. Boundary Yang-Baxter Equation. The variables g, /, g' ', /' = 1, 2, q satisfies 
the "admissibility conditions" g 7^ c, e; / 7^ 6, g; g' 7^ a, c; /' 7^ d, g'. Boundary conditions 
can also place further constraints on these variables. 

Fig. 6. Scattering processes described by the boundary unitarity condition (20); b 7^ 
a, c and c 7^ b. 

Fig. 7. Scattering processes for the cross-unitarity condition (22); b 7^ a, c and d ^ a,c. 

Fig. 8. Boundary bootstrap equation where states e and / must satisfy the "admissi- 
bility condition" and the boundary condition. 

Fig. 9. Physical poles of the boundary S-martix with conditions a 7^ b 7^ c in [Fig. (a)], 
and a 7^ 6, c in [Fig.(b)]. Here g£ is the boundary coupling constant for the particle A^. 

Fig. 10. Boundary S-martix elements for the free boundary condition with a, 6, and c 
all different. 
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